AD-768  208 


NON-LINEAR  WAVE  PROPAGATION 
R  c  S  .  Ri  v  lin 
Lehigh  University 


Prepared  for: 

Ballistic  Research  Laboratories 
September  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Fort  Royal  Road,  Springfield  Va.  22151 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  nay  be  obtained 
fro*  the  National  Technical  Information  Service, 

U.S.  Departaent  of  Conserce,  Springfield,  Virginia 
22151. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


<;  '--NR 


:'•*  utiiv  n.i  *«  ifi<  ni'*»i  f_Y/J  /  \js  &  IS  (7 

tMmuim  -  ■+*  +*4*4*—  tsom***ru*»*+  *%>»...  j iw^-iwm.  MB  <r<mjnw  iki  irrii  j»4  wcm^jwb.  twm+mmm  who—  i  iii  ■iri**lBfawMiwwiP<wlww>*,«  'Jwnwtiifcl»i>»t  i»in{iin» 

DOCUMENT  CONTROL  D^TA  *  R  &  D  " 

_ (Security  flat  *  Mention  ot  tltlo,  hotly  q  t  Abstract  and  Indexing  annotation  mu*t  be  entered  when  the  ovtratl  report  lx  ctnsaltlerh 

l.  ORiGlNA  riNC  ACTIVITY  (Corporate  vultiot)  REPORT  SECURITY  CLASSIFICATION  b 

« -i. :  _k  n.i _ I...  Unclassified  f 


Lehigh  University 
Bethlehem,  Pennsylvania 


J.  REPORT  TITLE 

Non-Linear  Wave  Propagation 


A.  DESCRIPTIVE  NOTES  (Type  ot  report  and  Inclusiv*  dates) 

Final  report 


S-  AU  TMOf»»ii  (first  no  me,  middle  Initial,  last  name) 

R.  S.  Rivlin 


6.  REPORT  DATE 

SEPTEMBER  1973 


9a  CONTRAC  T  OR  GRANT  NO. 

DAAD05-69nC-0053 

b.  PROJECT  NO. 


12b.  GROUP 


70.  NO.  OF  REFS 


ORIGINA-OII'S  RfiPORT  NUKfBERiS) 


BRL  CONTRACT  REPORT  MO.  120 


ob.  other  report  NO(S>  (Any  other  nunber*  that  may  be  aaalened 
thia  report) 


10.  DI3TRIDUTION  STATEMENT 


Approved  for  public  release;  d’Stribution  unlimited. 


II.  SUPPt.EMtaN'!  ARY  LOTUS 


!  13.  AU3TRACT 


I*.  SPONSORING  MILITARY  ACTIVITY 


US  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  MD  21005 


This  report  is  a  summary  of  the  work  accomplished  under  THEMIS  Project  No.  65, 
"Non-Linear  Wave  Propagation."  In  addition  it  lists  all  technical  reports,  and  the 
corresponding  papfers  published  in  the  course  of  the  project. 


Rsptodftf't  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  $  D*pOt!«»‘i*  of  Co*r>n**H» 
v  A  72151 


POM 

I  MOV  M 


1473 


M/4  fUrtiUKI  OO  F OH M  14/*.  1  J  AN  «4,  WHICH  IS 

/  /J  ODtOUKTA  row  ARMY  Uf*.  »  i 


UNCLASSIFIED 


H*«*urity  U>)'isi(icnUon 


Wave  Propagation 
Non  I  inear 
Elasticity 
l:luld  Mechanics 
Electromagnetic 
Wave  Propagation 
Shocks 

Nominear  Vibrations 


<r  {  ROUE  I  WT  ROLE  WT 


Security  Clun»lNcftUon 


/f£ 

CONTRACT  REPORT  NO.  120 

CENTER  FOR  THE  APPLICATION  OF  MATHEMATICS 

Lehigh  University 
Bethlehem,  Pennsylvania 

SEPTEMBER  1973 

FINAL  REPORT 
On 

(Project  THEMIS) 

Department  of  Defense  Contract  No.  DAAD05-69-C-0053 
"Non-linear  wave  propagation" 
monitored  by 

Ballistics  Research  Laboratory,  Aberdeen  Proving  Ground,  Md. 

prepared  by 
R.S.  Rivlin 


Approved  for  public  release;  distribution  unlimited. 


SUMMARY  OF  MAIN  ACHIEVEMENTS 


The  deformation  produced  in  non-linear  layered  media  when  sub¬ 
jected  to  dynamic  loads  has  been  analyzed.  Specifically  the 
concepts  of  non-linear  impedance  and  reflection  coefficients  of 
an  interface,  separating  different  non-linear  materials, has 
been  used  to  obtain  algorithms  describing  the  decay  of  a  pulse 
as  it  bounces  back  and  forth  in  a  slab  of  non-linear  material 
contained  between  two  other  materials.  Typically  the  deforma¬ 
tion  could  be  caused  by  the  arrival  of  a  shock  of  constant 
strength  in  one  of  the  outer  layers.  It  could  also  be  caused 
when  the  slab,  which  is  part  of  a  composite  material,  impacts 
on  some  other  non-linear  elastic  material.  This  study  is  poss¬ 
ible  because,  for  a  wide  range  of  model  materials,  whose 
constitutive  equations  approximate  those  of  man.y  real  materials, 
the  governing  equations  can  be  integrated. 

A  small  amplitude  finite  rate  theory  governing  resonant 
vibrations  in  crystals  or  in  gas-filled  tubes  has  been  developed. 
In  this  theory  the  effect  of  energy  radiated  from  the  face  of 
the  crystal  or  from  the  end  of  the  gas-filled  tube  has  been 
taken  into  account.  In  the  case  of  the  gas-filled  tube,  good- 
agreement  with  experimental  results  has  been  obtained. 

The  mode  of  propagation  of  a  long  gravity  wave  of  large 
amplitude  as  it  moves  into  a  region  where  the  fluid  is  sheared 
and  stratified  in  a  vertical  direction,  has  been  analyzed  in 
detail.  Such  waves  occur  in  the  atmosphere  and  tne  oceans. 

In  the  case  when  these  waves  o<_our  in  che  atmosphere  the  pre¬ 
dictions  of  the  theory  have  been  correlated  with  certain  aspects 
of  the  actual  behavior  of  severe  storms.  In  the  case  when  th*' 
waves  occur  in  the  ocean,  the  theory  has  been  used  to  model  the 
effects  of  strong  currents  on  tsunamis. 

The  techniques  used  in  classical  geometrical  acoustics 
have  been  extended  to  describe  the  propagation  of  large 


amplitude  waves  through  inhomogeneous  media.  This  theory  has 
been  used  to  describe  the  motion  of  a  plane  acoustic  pulse  in 
the  atmosphere/  where  the  ambient  stratification  is  due  to  the 
weight  of  the  air.  The  theory  has  also  been  used  to  describe 
the  propagation  of  a  tsunami  into  shallow  water. 

The  theory  for  the  propagation  of  plane  acoustic  pulses 
has  also  been  generalized  to  describe  the  propagation  of  non- 
planar  waves.  In  particular,  a  non-linear  ray  theory  for 
waves  propagating  in  elastic  materials  has  been  constructed. 

Theories  describing  the  propagation  of  large-amplitude 
short  pulses  in  viscoelastic  materials  have  been  developed. 

In  particular,  these  provide  techniques  for  the  study  of  the 
propagation  of  weak  shocks.  Also,  a  theory  has  been  developed 
which  describes  the  deformation  produced,  in  a  non-linear 
viscoelastic  material,  by  rapid  cyclic  loading  at  the  boundary. 

Work  has  been  conducted  on  the  existence  of  the  de¬ 
excitation  shocks  (6-shocks)  in  expanding  non-equilibrium 
nozzle  flows.  Situations  in  which  an  adiabatic  shock  (a-shock) 
is  embedded  within  the  6-shock  were  discussed.  It  was  shown 
that  the  de-excitation  process  cannot  lead  to  thermal  choking. 
Numerical  solutions  of  the  full  equations  were  also  obtained. 
These  solutions  confirmed  the  existence  of  the  6-sfcocks. 
Possible  asymptotic  limiting  states  were  deduced  for  a  general 
class  of  rate  equations. 

An  analysis  has  been  carried  out  of  the  validity  of  the 
sudden-freeze  approximation  for  supersonic  nozzl^  flows  of  a 
vibrationally-relaxing  gas,  when  the  initial  exv tation  level 
is  finite.  It  was  assumed,  in  agreement  with  e  ^erimentai  data, 
that  the  relaxation  time  changes  exponentially  with  the  trans¬ 
lational  temperature.  It  was  found  that,  in  this  limit,  the 
sudden-f reeze  method  does  provide  a  valid  first  approximation 
to  the  asymptotic  irozen  level.  It  was  shown  that  the  accu¬ 
racy  of  the  method  can  be  improved  by  introducing  a  modified 
freezing  point. 
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For  supersonic  nozzle  flows  of  a  condensible  vapor,  it  is 
often  observed  that  the  non-equilibrium  zone  downstream  of  the 
condensation  point  is  terminated  by  a  sudden  collapse  of  the 
supersaturated  state.  This  fact  has  led  to  the  concept  of  a 
condensation  shock.  Classical  shock  analyses  do  not  predict 
either  the  shock  position  or  the  thermodynamic  state 
immediately  downstream  of  the  shock,  nor  do  they  provide  infor¬ 
mation  on  the  conditions  under  which  such  shocks  will  occur. 

The  work  carried  out  in  this  area  attempted  to  answer  these 
questions  for  both  heterogeneous  and  homogeneous  nucleation. 
Conditions  for  the  existence  of  such  shocks  were  obtained  in  the 
high  activation  limit.  It  was  found  that  the  shock  is  pre¬ 
ceded  by  a  precursor  zone  whose  structure  is  dominated  by  drop¬ 
let  production.  A  simple  criterion  for  the  shock  position  was 
derived.  Further,  it  was  shown  that  the  asymptotic  state 
downstream  of  the  shock  corresponds  to  a  saturated  state  provided 
that  the  flew  does  not  choke.  Analysis  of  the  full  shock 
relations  implied  that  thermal  choking  will  not  occur  if  the 
flow  at  the  condensation  is  supersonic  with  respect  to  the 
frozen  sound  speed. 

Work  has  also  been  carried  out  on  the  propagation  of  plane 
electromagnetic  waves  in  a  general  isotropic,  non-linear,  non- 
dispersive  medium.  Formulae  have  been  derived  for  the  velocity 
of  propagation  of  the  shock  in  the  presence  of  constant  ambient 
electric  and  magnetic  fields.  It  is  found  that  the  propagation 
velocities  depend  on  the  electric  and  magnetic  fields  immed¬ 
iately  ahead  of  the  shock  and  on  their  magnitudes  immediately 
behind  the  shock. 

Formulae  have  also  been  derived  for  the  velocities  of 
propagation  of  first-order  electromagnetic  discontinuities.  It 
was  found  that  there  are,  in  general,  two  possible  propagat-Ion 
velocities  in  a  specified  forward  direction  and  correspondingly 
two  possible  polarization  directions.  Both  the  velocities  and 
the  directions  of  polarization  depend  on  the  magnetic  induction 
field  and  electric  displacement  field  immediately  ahead  of  the 
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discont.i  -uity.  Expressions  have  also  been  obtained,  in  the  case 
of  plane  'iscontinuities ,  for  the  manner  in  which  the  magnitude 
oi  the  discontinuity  changes  as  the  surface  of  discontinuity 
propagates . 

Simple  wave  theory  has  been  applied  to  the  determination  of 
the  manner  in  which  harmonics  of  various  orders  grow  as  an 
initially  sinusoidal  wave  propagates  in  a  non-linear  dielectric. 
This  method  is  more  powerful  than  those  usually  used,  in  that 
it  enables  us  to  calculate  the  amplitudes  and  phases  of  harmon¬ 
ics  of  all  orders. 

Our  results  agree  well  with  those  obtained  by  conventional 
methods  for  harmonics  of  low  orders. 

Work  has  been  carried  out  on  the  propagation  of  electro¬ 
magnetic  waves  in  a  non-linear  isotropic  material.  I u  has  been 
shown  that  small  departures  from  the  classical  linear  constitu¬ 
tive  relations  can  have  significant  cumulative  effects  on  the 
evolution  of  these  waves.  A  non-linear  far  field  theory  has 
been  developed.  Propagation  of  weak  electromagnetic  shocks  was 
also  discussed. 

Theories  have  been  developed  for  photoelastic  effects  in 
dielectrics  undergoing  finite  deformations  in  which  the  dielec¬ 
tric  constant  depevjds  on  the  instantaneous  values  of  the  defor¬ 
mation  gradients,  or  on  their  history.  The  latter  theory, 
which  may  be  expected  to  apply  to  viscoelastic  materials,  is 
developed  both  for  soiids  and  for  fluids. 
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